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Abstract 

In this paper the basic crack problem which is essen- 
tial for the study of subcritical crack propagation and 
fracture of layered structural materials is considered. 
Because of the apparent analytical difficulties, the pro- 
blem is idealized as one of plane strain or plane stress. 
An additional simplifying assumption is made by restric- 
ting the formulation of the problem to crack geometries 
and loding conditions which have a plane of symmetry per- 
pendicular to the interface. The general problem is for- 
mulated in terms of a coupled system of four integral 
equations. For each relevant crack configuration of prac- 
tical interest the singular behavior of the solution near 
and at the ends and points of intersection of the cracks 
is investigated and the related characteristic equations 
are obtained. The edge crack terminating at and crossing 
the interface, the T-shaped crack consisting of a broken 
layer and a delamination crack, the cross-shaped crack 
which consists of a delamination crack intersecting a 
crack which is perpendicular to the interface, and a 
delamination crack initiating from a stress-free boundary 
of the bonded layers are some of the practical crack geo- 
metries considered as examples. The formulation of the 
problem is given in Part I of the paper. Part II deals 
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with the solution of the integral equations and presen- 
tation of the results. 

1 . Introduction 

If one examines the evolution of a typical fracture 
failure in layered composites, one may invariably trace 
the initial cause to a localized imperfection which, from 
the viewpoint of loading, geometry and material strength 
corresponds to the "weakest link" in the medium. By far 
the most common forms of such imperfections which may 
have the potential of growing into a macroscopic crack 
and of causing eventual failure are the surface flaws, 
flaws in interlaminar bonding, and the lines of inter- 
section of the interfaces and free boundaries such as 
holes and other stress-free edges. Under cyclic loading 
and/or adverse environmental effects a surface flaw, for 
example, may grow into a part-through surface crack. 

Upon further application of the loads the surface crack 
may propagate subcritically through the entire^ m 
of the first layer. Following the path of least resis- 
tance, the crack may either propagate into the adjacent 
layer or grow along the interface. In analyzing the sub- 
critical growth of these delamination cracks as well as 
the cracks imbedded into homogeneous layers it is by now 
generally accepted that the stress intensity factor or 
some other parameter based >on the stress intensity fac- 
tor (such as the strain energy release rate in the case 
of delamination cracks) cart be used quite effectively as 
the primary correl ation parameter. Irt studying the frac- 
ture of layered materials the basic mechanics problem is 
then the calculation of the stress intensity factors 
along the crack front for all physical 1^^ exter- 

nal loads and crack geometries. 

The actual probl em i s a very compl i cated three- 
dimensional problem and at present seems to be analytically 
intractable. All the existing solutions are therefore 
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based on the two-dimensional or axisymmetric approxima- 
tions. Also to keep the analysis within manageable bounds 
in most of these solutions the medium is generally assumed 
to be infinite consisting of either semi -infinite spaces 
with or without a layer in between, or periodically stacked 
laminates. For example, the plane and axisymmetric pro- 
blems for a medium which consists of two or three differ- 
ent materials and which contains a crack parallel to or 
located at a bimaterial interface were considered in [1-4]. 
The problem of a crack perpendicular to the interfaces 
may be found in [5-9]. The problem of a T-shaped crack 
located on and perpendicular to the interface of two bonded 
half planes was discussed in [10]. The layered composite 
which consists of periodically arranged two dissimilar 
bonded layers with cracks perpendicular to the interfaces 
was considered in [11,12]. The effect of the elastic pro- 
perties and the thickness of the adhesive in bonded lay- 
ered materials was studied in [13]. The problem of an 
infinite medium which consists of periodic dissimilar 
orthotropic layers having cracks was studied in [14]. 

In this paper we consider a problem which is also 
idealized but at the same time is somewhat closer to the 
actual problem. It is a plane problem of two bonded iso- 
tropic infinite layers containing cracks of various orien- 
tations and sizes. The particular crack configurations 
which may be of considerable practical interest and which 
have been studied in this paper are shown in Figure 1. 
Unless one is dealing with a composite beam or a plate 
with through cracks, the idealization made in this paper 
for solving the problem is also rather severe not only 
because of the plane stress or plane strain assumption 

but also because of limiting the number of layers to two 

and assuming that the materials are isotropic. The com- 
posite laminates are, of course, multilayered and ortho- 
tropic. However, by a judicious choice of thickness 
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ratio and material constants, the solution given in this 
paper may provide valuable quantitative and qualitative 
information which may be useful in studying the real pro- 
blem. 

2. Solution of Differential Equations 

The fundamental problem which forms the basis of all 
the crack geometries shown in Figure 1 is described in 
Figure 2. Except for the interface cracks, the two iso- 
'tropic infinite elastic strips are bonded along they axis, 
the stnips contain arbitrarily oriented cracks along the 
X axis, and the problem is assumed to be symmetric with 
respect to the y=0 plane in geometry as well as applied 
Ibads. Furthermore, it is assumed that by a proper super- 
position the problem is reduced to a perturbation problem 
in which the crack surface tractions are the only external 
loads. Let the coordinate systems be selected as in Fig- 
ure 2 and let u^, V., (i=l,2) be the x and y coordinates 
of the displacement vector in the strips. The following 
differential equations must be solved for each strip under 
appropriate boundary and continuity conditions: 


(k-1 ) v2u+2( 
(k-1)v2v+2( 


9^U 

3X^ 


+ Aiv_) 

3xay ^ 


+ 

ax3y ay^ ' 


0 , 

0 , 


(la) 

(lb) 


y,here k=3-4v for plane strain and k=( 3 -v)/( 1 +v) for gener- 
alized plane stress, v being the Poisson's ratio. Because 
of symmetry , the problem will he considered for 0<y«>? only. 

Let the solution of (1) be expressed in terms of the 
following Fourier integrals: 

y(x,y)=| f [f(x,t)c6s y t+h(yit)Hnxt]dt, (2a) 

v(x,y)=^ [ [g(x,t)sinyt+k(y,t)c:os xt]dt. (2b) 



Substituting from (2) into (1 ) one would obtain a system 
of ordinary differential equations for the unknown func- 
tions f, h, g. and k. Solving these equations and obser- 


ving that u, V remain bounded as y-»-<» we find 
f(x,t) = -i- {[F-^ H + X t G]cosh xt 

+ 6 + X t H]sinh xt}, (3a) 

g(x,t)4 H + X t 6]sinh xt 

+ G + X t H] cosh xt } , (3b ) 

h{y,t) = (C+ytD)e"-y^ , (3^^ 

k(y,t) = (C+KD+ytD)e"-^^, (^^^ 


where C, D, E, F, G, and H are functions of the transform 
variable t and are unknown. Using the condition that c^^y 
must vanish for y=0 in both strips, which follows from 
the assumed symmetry, and defining A=tD, equations (3c) 
and (3d) may be replaced by 

h(y,t)=-^^ (^ - yt) (4a) 

yt)e-^*. (4b) 

From the stress-displacement relations and from (2), 
(3) and (4) the stress components in each semi-infinite 
strip may be expressed as follows: 


1 g (x,y) = -B- {[(F+xtG)sinh xt+(E+xtH)cosh xt]cos yt 
2]i XX 


f ( 1 -yt ) Ae'^^cos xt }dt. 


(5a) 
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-^yy(x,y)«|- J {£(F+2H+xtG)s1nh xt + (E+2G 

0 

+xtH)cosh xt]cosyt-(l+yt)Ae“'^*cos xt}dt, (5b) 

0CO 

axy^^»y^“l J {C(E+G+xtH)sinh xt + (F+H 

0 

+xtG)cosh xt]sin yt-yt Ae~^^sin xt}dt. (5c) 

Referring to Figure 2 we note that equations (2-5) are 
valid with u=iJi, k^icj, 0<y<~, 0<x=xi<2hi for strip 1 and 
with u-u 2 » k=k 2 » 0<x=X2<2h2 for strip 2. Thus, 

there are all together ten unknown functions, A^. , , F^. , 

G^. , and , (i=l,2) which may be determined from the fol- 


lowing boundary and continuity conditions: 

avxx( 2 hi,y)= 0 , a^j^y( 2 hi,y)= 0 , 0 <y<»iV ^ (6a, b) 

a2xx^0>y)“0» 0<y<», (7a, b) 

0<y<«> (8a ,b ) 

ui (0,y)=U2(2h2.y) , 0<y<a3, b 3 <y<», (9a) 

0ixx(o>y)=P3(y) > 33<y<b3 . (9b) 

Vi( 0 ,y) = V 2 ( 2 h 2 ,y) , 0<y<a3, b3<y<» , (Ida) 

aixy(0,y)=P4(y) . a3<y<b3 , (10b) 

V| (x^ , 0 )- 0 , 0 <x^<a^ , b^<X;j< 2 |i^fli;(i^l ,2) , (11a) 

0 jyy (x^ »0) = -p^ (x^ ) , a.<x^.<b^ , (i-1,2) . (11b) 


Substituting from (5) into the six homogeneous con- 
ditions (6)-(8), six of the ten unknown functions may be 
eliminated. The mixed boundary conditions (9)-(ll) would 
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then give a system of four dual integral equations to deter- 


mine the remaining four unknowns. The problem may also be 
reduced directly to a system of four singular integral 
equations by defining the following four new unknown func- 

tions : 


^-^v^(x.,0) = <l>^(x^) » (i®l»2) 

(12) 

-g|“Cvi(+0,y)-V2(2h2-0,y)]=«l>3(y) . 

(13) 

-:^Cux(+0,y)-U2(2h2-0,y)]=4)4(y)- 

(14) 


If we now replace the mixed conditions (9)-(ll) by (12)- 
(14) and observe that 

<l>.j(x^)=0, 0<x.<a^. , b|<X|<2h|, (i=l,2) , (15) 

<l> 3 (y)= 0 » <|> 4 (y)= 0 , 0 <y<a 3 , b 3 <y«» , (16a, b) 

from (2), (3), (4), (5)-(8), and (12)-(14) we obtain 

A|(t) = - j %^(s)sin ts ds, (i=l, 2 ), ( 17 ) 

"i 

and 

(Fi+2hit6i )sinh2hit+(Ei+2hitHi )cosh2hit=Ri (t) , 

(Ei+6i+2hitHi)sinh2hit+(Fx+Hi+2hitGi)cosh2hit=R2(t) , 

E2“Rs(1')» 

uiEi-U2C(F2+2h2tG2)sinh2h2t+(E2+2h2tH2)cosh2h2t]=R4(t) , 

yi(Ei+Hi)-y2C(E2+S2+2h2tH2)sinh2h2t 

+(F 2 +H 2 + 2 h 2 tG 2 )cosh 2 h 2 t]=R 5 (t) , 

Fi — 2 — H X - C (E2"^^'^^^^2'*‘2h2tH2 ) s i nh2h2t 

+ ( F2“^^^^^ H 2‘^2h 2fEI 2 ^ ^2h 2 i* ^“^6 ^ ^ V* 
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Ei+<i^ Gi-[(F 2+-^^ H2+2h2tG2)sinh2h2t 

+ (£2*^ *^^ ' ^” H 2 ) CO s h 2h 2 1 ] ® Ry ( t ) » (183-h) 

where Ri,..,R 7 are given in Appendix A in terms of 
<|) 4 , and the subscripts 1 and 2 in the unknown functions 
refer to the strips 1 and 2, respectively. 


3. The Integral Equations 

By solving the system of algebraic equations (18) 
and by using (17) it is seen that the stresses and dis- 
placements in the strips may be expressed in terms of the 
new unknown functions <|)i » . . »<t >4 only. The functions 
<|)4 may now be determined by using the boundary conditions 
(9b), (lOb), and (lib) which have not yet been satisfied. 
Thus, by substituting from (5b), (17) and (18) into (11b), 
after somewhat lengthy manipulatidiis we obtain 


V e ■ V 


s-x. s+x. 
®i ^ 


:^)c|,.(s)ds+ i .(x. ,s)(|..(s)ds 

^i V * 


=_ILU±J^)pi (Xi ) , i=l ,2, a^.<x.<b|, a4=a3, b4=b3,(19) 


where if a^?«0 and b^.?*2h^. the kernel s k . ^ , ( i -1 ,2 ; j =1 , . . ,4 ) 
are bounded functions in their respective closed domains 
of definition. The expressions for k^j are, of course, 
dependent on the solution of the System (t8)whic^ is quite 
cumbersbme and hencOv will not be given in this paper. The 
compl ete detail s may be found in D1 5] . 

Similarly, by substituting frorti (5a), (5c), (17) and 
(18) into (9b) and (10b), after ^e|IMting the dominent 
parts of the kernels we obtain 



, r-. : . , o:- 



4-r 

»»21 j=1 


f jk3j(y,s)+^ (s)ds.l^^P3i 




5 



( 20 ) 


a4“®3» b4*b3, 

■'83 

rb. 


ira2i 


1 Pk,j(y,s)+j(s)ds=|^^P4(y). 

J — I ^04 


where 


a 3 <y<t> 3 , 34=83, b4 = b 3 , (21) 


( 22 ) 


and the kernels ksj and k^|, (j=l,..,4) are given in [15]. 

From the definition of the functions <|>i , . . ,({.4 gi ven 
by (12)-(14) and the conditions (15) and (16) it is clear 
that for the imbedded cracks shown in Figure 2 the solu- 
tion of the integral equations (19)-(21) must satisfy the 
following singlevaluedness conditions: 

(s)ds=0, (i = l , . . ,4) , 34 = 33 , 64 = 63 . (23) 


4. Singular Behavior of the Solution 

For each crack configuration shown in Figure 1 the 
singular behavior of the solution of the integral equa- 
tions (19)-(21 ) and that of the stress state around the 
crack tips or the irregular points a^, b^., (i = l,2,3) may 
be examined by using the function theoretic method described 
in [16] (see also [17] and [18] for applications to crack 
problems). For ai>0, b 2 < 2 h 2 , and a 3>0 the dominant parts 
of the integral equations are uncoupled. Following [16], 
if we express the solution of (19) by 
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g .(,$) . . 

(b.(s)* » 94'<S<b., {j«l,2) 


(s-aj )'*^(bj-s) j 

0<Re(«j ,3j )<1 


( 24 ) 


substituting from (24) Into ( 19 ) , for the imbedded cracks 
it may easily be shown that [17,18] 




(25) 


In (24) g- is a Holder-continuous function and is nonzero 
at the ends a. and b., (j=l,2). Also, defining 

J J 

i|)(s) = (!)4(s )+ i <|)3(s) > 

the integral equations (20) and (21) may be combined as 


1 

iri 


''3 ^ 

as 


If we now let 

1J)(S) = 


ds-Yi|'(y) = a^”CP4(y)-i Paly) 3 
+ bounded terms, a 3 <y<b 3 . 

93(s) 


(27) 


(s-a3)“Mb3-s)®5 


, 33<s<b3, 0<R(a3 ,&3 )<1 J (28) 


from (27) following the precedure outlined, for example, 
i n [1 8] we obtai n 

a3=^-io), 3 3=^1“» <^“2^ ’ (29) 

where 93 is again a H61der continuous function and is non- 
zero at the ends a 3 and b 3 * 

For the limiting cases of cracks intersecting each 
other, free boundaries, or the interface the kernels k^j 
which appear in (19)-(21) are no longer bounded for all 
values of their arguments. In such cases to determine the 
correct singular behavior of the functi ohs i>i » ♦ * » 

through a proper asymptotic analysis these unbounded parts 
of the kernels must be separated and must be taken into 
consideration in the application of the function theoretic 
methods to the integral equations. The typical crack tip 
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behaviors which do not conform to the standard singulari- 
ties expressed by (24), (25) and (28), (29) are discussed 
below. 

(a) Edge crack, bi*2hi. 

Consider, for example, the problem described by 
Figure 2 in which bi=2hi, ai>0, b2<2h2» a2>0» a3>0, b 3 <«>. 

In this case of all the Fredholm kernels appearing in (19)- 
(21) only kii(xi,s) exhibits singular behavior. Analyzing 
the integrand which defines kn we find that 

-00 00 

kii(xi,s) = J CKii(xi,s,t)-Kn(xi,s,t)]dt 

^ 00 00 b S 

+ ] Kii(xi ,s,t)dt=kii+kii , (30) 

* 0 

kii(xi»s)=| Kii(xi 

^ o 


6(2hTXi ) _ 4(2hi -xi 

(Xi+s-4hi)^ (xi+s-4hi)^ 


(31) 


where K^i is the asymptotic value of Kn for t->~, xi-^2hi, 
and s-»-2hx. Now substituting from (30), (31) and (24) 
into (19) and using the technique described in [18] we 
obtain the following characteristic equations giving the 
powers of singularity aj and Bi* 

COtirai = 0, (32a) 


gi(2hi) 

sinirBi 


[2(1 -Bi ) 2-1 -coSTr6i3=0. 


(32b) 

\ 


The acceptable roots of (32) are 
oj — 0.5 ,Bx®0, 

which are the known results [19]. Powers of singularity 
identical to (33) are obtained for a2=0, b2<2h2- 
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(b) Crack tip terminating at the interface, ai=0. 
In Figure 2 let ai=0. In this case too ku is the 
only Fredholm kernel which exhibits singular behavior. 
Let kji again be expressed by (30). The asymptotic part 
of the integrand may be separated as follows: 


where m=ui/y2* Noting that for ai=0 in (19) the term 
l/(s+Xi) also becomes unbounded for s->-0, Xx-^0, the addi- 
tional singular part of the kernel may now be expressed as 


k ! 1 ( X 1 , s ) I Kt 1 ( X 1 , s , t ) d t 

0 2 
Cii ^ ^12X1 ^ Cx 3 Xx 

S+Xx ^(s+Xx ) ^ ^ (s+Xx ) ^ ’ 


_1 1+Kl 

" 2(l+mic2) 


3(m-1) 

2 (m+K X ) ’ 


m=v x/P2 


P - 6(n»-1 ) 

^12 m+Ki 


Cl 3“ 


4(m-l ) 
m+K X 


(35a) 


(35b) 


The kernel l/(s-x) combined with (35) constitutes a gener- 
alized Cauchy kernel. Following the procedure of [18], 
from (19), (24), (30), and (35) the characteristic equations 
which determine ctx and Bx for ax=0 may be obtained as 

coti:Bx"0» (36a) 

COSirax + Cxx Cx2“l ^xs — I ^ ^ ^ = 0. (36b) 

which are the known results [5]. 

Similarly, for b2=2h2» a2>0, ax>0, bx<2hxj and a3>0 
we obtai n 
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k|2(X2 »S ) = 


C 22 ■ C 22 C2h2'*-X2 ) ■ '^2 S C2h2'"X2 ) ^ 


4h2"^2"S (.4h2-X2*'S|^ (4h2-X2-S2 ) ^ 


1 ^ / n m(T+IC2) 

so^iiir s'(i+Vi'r ■ 

P . 6(m-1) . , 4(1-ni) 

T+m^r ^^2 3 i+n,K2 * 


( 37 ) 


From (37), (24), (30) and (19) the characteristic equations 
are found to be 

COtira2*0 , (38a) 

CO Sir 02 “ ^21 “ 32^22 ” 2 ( "^ 2 ) ^2 3 (38b) 

which are again the known results [5,6]. 

(c) Crack intersecting the interface, ai=0, b2=2h2- 
Referring to Figure 2, if ai=0, b2=2h2» a3>0, 
bi<2hi, and a2>0, in (19) not only kn and k 22 but also 
the coupling kernels ki 2 and k 2 i exhibit singular behavior. 

The singular parts of ku and k 22 given by (35) and 
(37). Using a notation similar to (30), for the singular 
parts of the coupling kernels we obtain 

f. c + f h rl+ct ?h tne‘‘'^^2h2-s+XiJ , 

Ki2(xi .s ,t)-[jj^:j:^(xit-2)“i+,nK2(^st-2h2t)]e ' ffM 


,.s _A_ dll ^ di2Xi 

kl 2 (Xl.s)- Xi-sizhT * (x,-S + 2 h 2)2 


I - 3(l+'^i) j 

■ rtm+TiT’ 


-I+lL _ 

^2“ni+K;j “ l+mK2* 


(39) 


Kji(xx.s.t)-[|f^{l-2st) 

. |i^(|+x,t-2hat)]e-‘(2''4-X4*^^ 

l<2l(X2.s)“ }h2 + S-X2 ' (ihi+S-Xjl' ’ 
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-Mil+lL) ELL+liy 
2 ^ "2 (1 +IUK2I *'2(jn+Ki) * 


d 


_ ntn f *^2 ) ' ih(T-^'^2 ) 

ifl^ici T +WK 2 


(40) 


Substituting now from (24) into (19), noting that at the 
irregular point Xi=0, X2®2h2» «i“32> using (35), (37), 

(39), and (40) we find the following Gharacteristic equa- 
tions to determine ai*B2» > ^nd 02: 


GOtTrei=0, GOtira2=0, 

[GpSiroi+Gn + aiGi2‘*'^“l ( 1 1 ) C 1 3 ] I- COSirai -O2 1 "Cti G22 

-^1 (l+ai)G23]+(dii+aidi2)(42l'^«l‘^22)"0» (41a-G) 

whiGh are the known results [6]. To obtain (41) it is 
assMmed that gi(bi)?«0, 92( Equation ( 41 g) is the 

expression of vanishing determinant of the linear homo- 
geneous algebraiG system in 9i(0) and 92(2h2)« This indi- 
catss that the constants 9i (0) and 92 ( 2h2 ) 3re not inde- 
pendent and are rel ated by 

Ccosivai + G 1 i + aiG 1 2‘*''^ 1 ( ^ ) G 1 3 3 ( 0 ) 

( d 1 1 +ai d 1 2 ) 92 ( 2h 2 ) “0 . (42 ) 

Condition (42) replaGes ons of the si ngl e-yal uodnoss con- 
ditions (23) in solving the system of singula'^ integral 
equations. 

For a3>0 the results given by (29), (32), (36) , (38) , 
and (41) cover all crack configurations shown in Figures 
(la) to ( 1 h ) . 

(d) T-shaped or cross -shaped oraoM- 

The problem becomes somewhat morocompli Gated if 
the interfaGe craGks interseGt the GraGhs WhiGh are per- 
pendiGular tp the boundaries, Oonslder., for example, the 
cross-shaped Graok shown in Figure 1(d) for which a^^O, 
hi<2hi, a2>Q, b2=2h2» 33=0, and |3<^ (see Figure In 

this Gase the point of interSeGtion (xi=50, yep) is an 









irregular point common to all four integral equations given 

by (19)-(21). Examining the asymptotic behavior of the 

kernels it is found that >2 >3 ,4) , and k 3 ^» 

k (i*l ,2) become unbounded as the variables Xj , X 2 , 

and y, approach the common irregular point together with 

s in pairs. Singular part of each one of these kernels 

may again be separated in a straightforward manner by 

analyzing the asymptotic behavior of the related integrands. 

After separating the singular parts k.j of the kernels, 

the integral equations (19)-(21) may be expressed as 

2h2 

+ k! 2 (xi.s)+ 2 (s)ds 


+s I kjj (xi ,s)<^j (s)ds+Pi{xi )=- ^ » 

0<Xi<bi, pu 


k^ (X2,s)(|)i(s)ds+| k|2{x2.s)](|)2(s)ds 


D3 

+ Z j k|j(X 2 »s)<|)j(s)dS+P 2 (X 2 )=- 


ir(l+'^2 


:(X2), 


a2<X2<2h2 


_2 

ira2i 


'^*k|i(y.s)+i(s)ds-,5|^ I k|2(y,s)*2(s)ds 

a 0 


b 3 

t 7 I (7i3r^i77)K(s)<is+Y't'3(y)*'’3(y) 


1+^1 

321V1 


P3(y)» 0<y<b3, 


b f ^ 

I *k4i(y.s)fi(s)ds+:j^ I kj2(y,s)4.2(s)ds 
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b3 


1+»^i 

aaiwi 


P4(y)» 0<y<b3, 


( 45 a-d) 


where the functions P,^, (k=1 ,4) represent all the remain' 

ing terms in the integral equations corresponding to the 
bounded kernels, kfi, ki2» ^22 given by equa- 

tions ( 35 ), ( 39 ), ( 40 ), and (37) , respecti vely , and 

AiXi AaXitxi + 5s^) 

kf3(xi,s) = ■— V - 

s2+xi (s2 +x2)2 

g AxS A2s(3s2.x2) 

ki4 (Xi ,S ) = - -- “ / 2 , 2\2 ’ 

S^+X^ (s2+X^)^ 

g ^ ^ A3(2h2-X2) A4(2h2 “X2)[(2h2“X2)^+5s2] 

k 23 (X 2 .s)-^^^^^^^_^^^^ - (-s2+(2h2-X2)2]2 

A3S A4SC3s2-(2h2-X2)^] 

k24(X2.s)— s 2 +( 2 h 2 -X 2 )^ ^ Cs2+(2h2-X2)n2 

2 AiS 2 A 2 s( 3 y 2 -s 2 ) 
k3i(y.s) = - — — -+ , _ ox, " ’ 

s2+y2 (s2+y2)2 


kf 4 (^1 )=- 


k|3(X2,s)=- 


k24(X2>^)~'“ 


k|i (y .s) = - 


k 1 2 ( y > s ) = 


2A2(2h2-s) 2Ai(2h2-s)C3y^-(2h2-s)2] 

'y^+{2h2-s)2~ [(2h2-s)2+y2]2 


2Aiy 2A2y(y2-3s2) 

ktfi (y >s)“ ' “ ’ 

S^+y2 (s2+y2) 2 


^ y^+(2h2-s)^ C (2h2-s )2+y2]2 
^1" 2 ( 1 +mic o') * 2{m+K i ’ ^ 2(m+Ki ) 


2A2y 


(s2+y2)2 

2Aiy[y2-3(2h2*s)23 


(46a-h ) 
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(47) 


A. ™ • 

Note that the T-shaped crack shown In Figure 1(i) is 
a special case of the cross-shaped crack for which bi=0. 

In order to simplify the manipulations in the asymp- 
totic analysis we first define 

x*=2h2’“X2» s =2h2"S, b =2h2“a2> ^2 )*'<l>2(s)» (48) 

in the interval a2l(x2»s)<2h2 and then, for convenience, 
drop the superscript (*). Thus the origin of the coordinate 
systems becomes the common irregular point of the integral 
equations (45), Noting that at s=0 the unknown functions 
<1)1,.., (|)4 must all have the same singular behavior, we let 

J 

gj(0)/0, (j®l>**>4)* (49) 


Using (49) and the procedure outlined in [16,17,18] one 
can establish the following asymptotic relations for the 
relevant singular integrals in the close neighborhood of 
the end point x=0, y=0; 


rb 


<I>(S) j.. g(0)cotir_g 

s-x bPx« 


irj S + X b^x“sinira 

0 




1 f^x<t(s)\i„ aq(0) Lf a(a+l)g(Ql... . 

7 ] ^ bVsli,;T’ 2bex“s1n.a 


1 


s<fe(s) g(0) . If g(Q.)— 

-2b®x«sin^’ ^JoS2+y2 2bVcos^ 


^i^(s 2 +x 2)2 "4b^x“sin^ 

- 17 - 


f'sxiiiLl A-- qg(o) - l\ nhhl 

(50a-h) 



Substituting now from C50) into (45) and using (46) 
and (48), tKe leading terms of the integral equations around 
the end point a*Q, y«Q ma/ be expressed as 

9l(0) t-, \ 

— 5 ^ ~£c0Siro+Cai+oCi2'’'Ci 3“vl+a;/Z j 

bi^^x^si nira 


+ 


9z(0) 


93(0) 


2b3®*x“ 


TTa 


[ei i+ei2(l -a)/2] 


XjCOS^ 


94(0) 


OK 64 „<x_,* _ira 

Zb 3 ^XiSin-^ 


(f 1 i+f 1 2<*/2 ) ”Fi (xx ) ^ 


gx(0) 92(0) 

(d2i + qd22)~^ ft— C-C0Sira+C2i + aC22 

bi n-ira b2^^X2Sin7ra 


+ G2 3«(l‘‘’a)/2]+- — — rrCe2i ■*■ 022(1 -a)/2] 
Zba^Sx^cos^ 


94 ( 0 ) 


OK S 4 „a_j_Tra 

Zb3 ^X2Si n~Y 


[f21+f22»/2]=F2(X2) > 


2gi(0) 


L S 1 a . IT a 
321^1 sin-2 


[Ai+(l-2a)A2]- 


292(0) 


^21^2 S 1 


c A 2 + (1 -2a) All 


T93(0) 94 ( 0 ) 

+ : — + 


bs^^ya b3^'+y®simra 


(1 +C0S7ra) = F3 (y ) » 


29i(0) 


a2it>i ^^y“cos-^ 

93 ( 0 ) 


CAi-(l-2a)A2]+ 


292 ( 0 ) 


32 192 COS-^ 


CA2-(l-2a)Ai] 


b3^ 3 y®siniro 


Y94 (0 ) 

-(cosira-l ) - , Q • - ~ F 4 (y) » 
t>3®V 


(51 a -d) 


where the constants c^j and d^.j are 9 iVqn by (35b), (37) 
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(39) , and (40), 

eiji*(Ai+A2) » ei2?*4A2 > e2i"’‘(A3+A4) , e22**'4A4 > 
^11“3A2-Ai, fi2*’‘^Ai , f21**'(A3+A4 ) , f22*^A4, (52) 

and the functions Fi,..,F 4 contain all remaining terms which 
are bounded at the end point x=0, y=0. If vfe now multiply 
both sides of equations (51a)-(51d) by xi» ^ 2 * y”’ y » 

respectively, and let Xi->-0, x 2 ^ 0 , and y-»-0 we obtain a system 
of algebraic equations for gi (0) , . . ,g 4 ( 0 ) of the following 
form: 

E B,/g.(0) = 0, k=l,..,4 (53) 

j = l ^ ^ 

where the coefficients are given by (51). Since gji (O) . 
g4(0) are assumed to be generally nonzero, from (53) the 
characteristic equation which accounts for a is obtained 
to be 

det|B^j| = 0, (k,j) = (l ,2,3.4). (54) 

After some simple manipulations from (54) it is found that 

-^Ccosira+l-2(a-l)2] = 0. (55) 

SI riTTa^ 

which is identical to the result found for the edge crack 
(see equation 32b). 

For the T-shaped crack shown in Figure l(i) = 
the problem is formulated by the integral equations (45b)- 
(45d), and the characteristic equation for a may be obtained 
from (51b)-(51d) as follows: 

det|B,j|-0. (1.j)-(3,3.4). (56) 

It can be shown that ( 56 ) too reduces to (55). 

The foregoing asymptotic analysis is restricted to 
the singular behavior of the solution at the common 
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Irregular point x«Q, y»Q, The analysis given in the pre- 
vious sections indicate that and the characteristic 

equations for gi,&2»&3» and « are uncoupled. Hence givgz, 
and 33 are given by (25) and (29). 

One should note that theoretically the solution of 
the system of singular integral equations (19)-(21 ) con- 
tains four arbitrary real constants [16]. In the case of 
nonintersecting cracks these constants are determined from 
four singlevaluedness conditions of displacements given by 
(24). However, in the case of intersecting cracks, such 
as T or cross-shaped cracks considered in this section, 
kinematically it is clear that there is only one single- 
valuedness condition, namely (see the definitions (12) and 


(13)) 

2I12 

* 

<p2(xz)<i>^z 
a 2 


ba 

<j>3 (y)dy+ 


fb 1 

<|)i (xi )dxi = 0. 



(57) 


The additional three conditions which are necessary for a 


unique solution of the system of integral equations are 
provided by (53). Note that with (54) satisfied, (53) 
gives three equations relating the end values g x(0 ),.. , 
g4(0 ). In the case of T-shaped crack 4>i = 0, (53) consists 
of three homogeneous equations and gives two independent 
conditions relating g2» 93 and 94. 


5. Stress Intensity Factort 

A careful exami nation of the i ntegral equations (19)- 
(21) would indicate that at a given irregular point if the 
displacement derivatives have a singularity of power a, 
then the stress state is also singular having the same 
power a. In applications it is important to know not only 
the power of stress singularity but also the coefficient 
of the singular term in the plane Of projected (or conjec- 
tured) crack extension. This coefficient is known as the 
strength of the stress singularity or the stress intensity 
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factor. TKese coefftctents can be evaluated in terms of 
the dtspl aceiijent dert^vatt^es . . ,4>»t ti/ simply observing 
that in the integral equations 09) •‘C21) the e^cpressions 
on the left hand side represent the related stress compo- 
nent outside as vfell as inside the cuts Ca.j» b^), Ci 
(a) Imbedded Crack. 

Consider the crack in material 2 and let 0<a2<t>2<2h2 
(Figure 2). Defining the sectionally holomorphic function 

b2 




a2 


^i(s) 

s-z 


ds , 


(58) 


and using (24) and (25) it can be shown that 

92(^2)®’^^^^ T 92(1^2) ^ 

(b2-a2)*^ (z-a2)^ (b2-a2)‘^ (z-b2)* 


F(z) = 


+ 0((z-c)“), (c=a 2 or b 2 » Re(to)>-%). 


(59) 


In (58) since F(z) is holomorphic outside the cut, from 
(58) and (59) it follows that 

‘’ %2(s) jS , 1 

S-X2 ^ 27 (b2-a2)2 


a2 


X2 ' - VD2-«2^‘ (32-X2)' 

^ 2 (^ 21 , L_ 3 _ + 0 (|x 2 -cr), (c=a 2 or b 2 ). 

(b2-a2) * (X2-b2) ^ ^gQy 


Referring to (19) and observing that the equation gives the 
expression for a 2 yy(x 2 , 0) for 0 <X 2 < 2 h 2 and that the terms 
containing the kernels which are bounded at X 2 =a 2 and X 2 =b 2 
would have no contribution to the singular behavior of the 
stress state at a 2 and b 2 » the asymptotic behavior of o^yy 
may now be obtained from (60). Thus, defining the stress 
intensity factors by 

k(a 2 )“ lim C2(a2-X2)3^ ^ 2 yy^^ 2 . 0 ). 

X2->-a2 



kCbz)* ItJD I2,Cx2^bi2D^ ^2yy^'^2»Q)» 

X2-»-b2 /V 


C61a,b) 


vte find 


4ii 2 , SaCszI . ^ ■^^rr^r \-i^ / \ 

T 4>2(^2) » 

^ICb2-a2)/23^ XjW 


4V2 

k(b2)“ ■* T+r;^ 


.a^Cba) 


2 E(b2-a2)/23' 


4 VI 2 1 

»- 1 im ■jV-"E 2(b2"X2)3 ^ <^ 2 (^ 2 ) 
X2"^^2 ^ 


(62a, b) 


(b) Crack Terminating at the Interface. 

If the crack tip touches the interface the stress com- 
ponent of primary interest is the "cleavage" stress in the 
adjacent medium. For example, let a2®0, b2=2h2> ai>0, and 
33>0 (Figure 2). In this case the first equation of (19) 
may be wri tten as 

2h2 

a^yy(xi .0)4 I 

0 ‘ ' 


0<Xi<2hx , 


(63) 


where <}i 2 is given by (24) with 02 defined by (38b), kx 2 is 
the singular part of ki 2 given by (39), and Hi(xx) 

essentially represents all the remaining bounded terms. 
Defining again the sectionaTly holomorphic function F(z) 
by (58) and substituting from (39) and (24), for the lead- 
ing terms around xi=0 we find 


1 

7T 


I** % 2 fs)ds : 92 1 
xi-s+2h2' (2h2)\lnxS2 
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0 



Oij '«Su+2tL2 i 


6292(^^21 ' I " 
C2fl2l^5tnTr 82 


C64a,b) 


If we define tfie stress Intensity factor by 

kCZbz) = <^ivvC5(i,0), 

xi-^0 

from (63), (39) and (64) we obtain 

4yi d22'*'32^12, 92(2^2) 
k(2h2) = 1+^^ simr82 /H 7 ~~ 


(65) 


(66) 


It should be noted that the stress components in the small 
neighborhood of the point xi=0, y=0 may be expressed as 

k(2h2) , 

*^kij(r»®)' -— 37 ^kij^®^’ O<0<ir. i,j = r,0), 

'^2 ■■ ' (67) 

where the functions f,^^j are given in [5]. 

(c) Crack Crossing the Interface. 

In the point of intersection of the crack and the 
interface even though one may again define a single stress 
intensity factor, from the viewpoint of applications it is 
more convenient to express the asymptotic forms of the 
normal and shear stresses along the interface separately. 
For example, let b2=2h2» ai=0, a2=0, bi<2hi, and a3>0 
(Figure 2). In this case the interface stresses are given 
by (20) and (21 ) , whi ch , for small values of y, may be 
expressed as 

2yj 

Zhz 

+ f kfaCy »s)«|»2(s)ds] + H2(y)» 
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2>1 

■r^yCQ.y) > j^fc„0y,s)4Us)ds 

Zhz 

+ I ^ 2 Cy'»^U 2 Csy<is 3 + HaCy), (68a, b) 

0 

Hhere the singular fcernels are given fay ( 46 ) and the func- 
tions H2 and H3 represent the reTnaining bounded terms. Sub- 
stituting from ( 24 ) and ( 46 ) into ( 58 ) and a similar Cauchy 
integral with the density ^2. and using the results given 
by ( 50 ), from (68) one may easily obtain the asymptotic 
expressions for and o^^,. The stress intensity factors 
for o and a may then be defined and evaluated as follows: 

XX 

k„„(0) = lim y“ < 7 xv( 0 »y) 
y->-0 

2ui [Ai + (l -2a)A2]cji(0) CA2+(T-2a)Ai]g2{0) , 

= ___L_;_ — — + — — — — }, 

/bi sin(iTa/2) / 2 FT 2 sin(TTa/2) 

k^^(O) = lim y“a (0,y) 

xy y-»-0 

2 ui [Ai-(l-2a)A2]gi(0) ^ [A2-(l-2a)Ai]g2(Q) ^ 

^ /F^ COS(Tra/2) / 2 F 2 C0S(TTa/2) 

( 69 a, b) 

where 0=01=62 is obtained from the characteristic equation 

( 41 c) and the constants Ai and A2 are given by ( 47 ). 

(d) Interface Crack. 

In the case of an interface crack substituting from 

(28) into (27) it can be shown that the asymptotic behavior 

of the contact stresses in the small neighborhood of the 

crack tips is of the following form C2]: 

'^ixy ” w'^fy ) * 
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( 70 ) 


63 «3 

WjjCy) * (y^ba) Ca3*^y) . » y<9 3» b3<y, 

where M Is a bounded functton. Since 03 and 53 are cotnplex, 
one may define the stress Intensity factors as follows: 

kiCb3) + 1k2(b3)*ltin ,y) + 1 (0 ,y)J , 

y-^hs 


ki(a3)+ik2(a3)*liin w (y)Ecr, (0,y) + ia^j^ (0,y)]. 

y->-a 3 

{7Ta,b) 

In terras of the density function i{'(y) defined by (26) and 
(28) these complex stress intensity factors may be expressed 
as E2] 

ki(b3)+i k2(b3) = -g3(b3)T+7^ ’'1 - . 

ki(a3) + 1k2(a3)=g3(a3)Y^:^ (72a, b) 

where gs(y) is obtained from the solution of the system 
of integral eguations (19)~(21). Note that because of the 
definition (71) the dimension of the stress intensity fac- 
tors is az rather than the conventional a/T . 

After calculating the stress intensity factors the 
strain energy release rate for the crack propagating along 
the interface may be obtained from 12] 


3U _ ird+'^i; 
3 ^ 3 " 2yia2it3 


(k2 + kb 


1 3= (b 3 -a 3 ) /2 . 


(73) 


References 

1. F. Erdogan and G.D. Gupta, "The stress analysis of multi 
layered composities with a flaw", Int. J. Solids Struc- 
tures , Vol . 7 , pp . 39-61 , 1 971 . 


- 25 - 



2 . 


3. 


4. 


5. 


6 . 


7. 

8 . 

9. 


10 . 


n 


12 . 


13. 


14 . 


F. Erdogan and S.D. Gupta, "Layered composites with an 
Interface flaw", Int. J. Solids Structures, Vol. 7, 
pp . 1 089 •«^^ 1 Q7 , 1 971 . 

K. Arin and f. Erdogan, "Penny-shaped crack In an elas- 
tic layer bonded to dissimilar half spaces , Int. d. 
Engng. Science, Vol. 9, pp. 21 3-232 , 1971 . 

F. Erdogan and K. Arin, "Penny-shaped interface crack 
between an elastic layer and a half space , Int. J. 

Engng. Science, Vol. 10, pp. 115-125, 1972. 

T.S. Cook and F. Erdogan, "Stresses in bonded materials 
with a crack perpendicular to the interface , Int. J. 
Engng. Science, Vol. 10, pp. 677-697, 1972. 

F. Erdogan and V. Biricikoglu, "Two bonded half planes 
with a crack going through the interface", Int, J, Engrg. 
Science, Vol. 11, pp. 745-766, 1973, 


D.B. Bogy, "The plane el 
metrically loaded crack 
Engng. Science, Vol. 11, 


astostatic solution for a sym- 
in a strip composite", Int. J. 
pp. 985-996, 1973. 


G.D. 

Int. 


Gupta, "A layered composite with 
J. Solids Structures, Vol. 9, pp. 


a broken laminate" 
1141-1154, 1973. 


J 


K. Arin, "A note on the fracture of laminated composites , 
Letters in Applied and Engineering Sciences, Vol. 3, 
pp. 81-85, 1975. 

J.G. Gpree and W.A. Venezia, "Bonded elastic half planes 
with an interface crack and a perpendicular intersec- 
ting crack that extends into the adjacent material , 

Int. J. Engng. Science, Vol. 15, pp. 1-27, 1977. 


F. Erdogan and M. Bakioglu, "Fracture of plates which 
consists of periodic dissimilar strips'^ Int. J. of 
Fracture, Vol. 12, pp. 71-84, 1976. 


F. Erdogan and M. Bakioglu, 
in layered materials", Int. 
pp. 739-749, 1977. 


"Stress-free end problem 
J . of Fracture, Vol. 13, 


M.R. Gecit and F. Erdogan, "The effect of adhesive 
layers on the fracture of laminated structures ,J 
Engineering Materials and Technology^ ASME, 
100, pp. 2-9, 1978. 


of 
Vol . 


F. Delale and 
with cracks", 
15, 1979. 


F. Erdogan, "Bonded orthotropic strips 

to appear in Int. J. of Fracture, Vol . 


- 26 - 



15, Mlng-Che Lu, "A Gojopostte of Two Bonded Strips Con- 
taining Perpendicular Cracks and an Interface Crack", 
Ph.D, Dissertation, Lefilgti University, October 1978. 


16. N. I. Muskhel Istivll 1 . WguTaf Integrar ^ 

P, Noordfioff, 6inonlngen , The Netfterl ands , 1 953. 

17. F. Erdogan, "Complex function Technique", Continuum 
Physics, Vol . II , A.C. Erlngen (ed), Academic Press, 
pp. 523-603, 1975. 

18. F. Erdogan, "Mixed Boundary-Value Problems in Mechanics " , 
Mechanics Today , Vol. 4, S. Nemat-Nasser (ed.), Pergamon 
Press, Oxford, pp. 1-85, 1978. 

19. G.D. Gupta and F. Erdogan, "The problem of edge cracks 
in an infinite strip”, J. Applied Mechanics , Vol . 41 , 
Trans. ASME,pp. 1001-1006, 1974. 


4 





- 27 - 




Appendix A 

Th.e functions Ct 1 , . . »R 7 Cfl ♦ 

R2(t)= i[l-(2hi+s)t]e‘*('""'^^' 

1 ^ai 

-[l-C2hi-s)t]e’^^^^^“^he^^^^'|)i(s)ds, 


R3(t)= 


1 


2tse"^^(|)2(s )ds , 


R4(t)= ^^[^^2tse"^%i(s)ds-;^| %C(2h2+s)e"^^'^'^2+s ^ 

K 1 ^ I J Z J 

3l 

-(2h2-s)e“^^^^2"^he^^2%2(s)ds. 

R5(t) = — %f^^{n-(2h2+s)t]e"^^'^^2+s) 

•'32 

-Cr-(2h2-s)t]e‘''^(^^2-s )}e2li2t^2(s)ds , 

R6(t)-^f C2h2+s)t]e-*'2*'2**> 


2 ' “a 


+ ( 2 h 2 -s)t]e*^^^^ 2 "^H<l> 2 (s)ds + f H(s)sints ds, 
^ ^ . *^3 3 

R7(t)=;j^| '[2ts-(Ka+l)3e'*%i(s)ds 

. 1 f'’^{[i|±L . C2h^ts)tje-‘(2R2+M 

<2+ ' Ja2 

,£l|ii . (2hj-s)t]e-‘^2ha-s)j^^(3) „3 

ft>3 

+ (j )3 (s )cos ts ds . 

J So 
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Figure 1 


Crack ‘configurations considered in the paper. 










Figure 2. 
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